Mon. Not. R. Astron. Soc. OOP. fUfMl (2012) Printed 18 September 2012 (MN WT^ style file v2.2) 



Hot Gas, Cold Gas and Sub-Halos in a Lyman-o^ Blob at 
Redshift 2.38 

"(Ni" Paul. J. Francis^'*', Michael A. Dopita^'^, James W. Colbert^, Povilas Palunas^, 

^ Claudia Scarlata^, Harry Teplitz^, Gerard M. Williger^''''^ and Bruce E. Woodgate^ 

^Research School of Astronomy and Astrophysics, The Australian National University, Canberra ACT 0200, Australia 
'^Astronomy Department, King Abdulaziz University, P.O. Box 80203, Jeddah, Saudi Arabia 

Spitzer Science Center, California Institute of Technology, Pasadena, CA 91125, USA 

Las Campanas Observatory, La Serena, Chile 

Minnesota Institute for Astrophysics, School of Physics and Astronomy, University of Minnesota, Minneapolis, Minnesota 55455, USA. 
^Department of Physics & Astronomy, University of Louisville, Louisville, KY 40292, USA 
^Lab. Lagrange, U. de Nice, UMR 7293, 06108 Nice Cedex 2, France 

^Inst. for Astrophysics and Computational Sciences, Catholic U. of America, Washington DC 20064, USA 
^NASA Goddard Space Flight Center, Greenbelt, MD 20111, USA 



Draft 6 Sept 2012 



o 

u 

Q ! ABSTRACT 

\^ ■ We present integral field spectroscopy of a Lyman-a blob at redshift 2.38, with a 

^ \ spectral resolution three times better than previous published work. As with previous 

^ . observations, the blob has a chaotic velocity structure, much of which breaks up into 

multiple components. Our spectroscopy shows, however, that some of these multiple 
components are extremely narrow: they have velocity widths of less than 100 km s~^. 
Combining these new data with previous observations, we argue that this Lyman- 

■ a blob resides in a dark-matter halo of around lO^^M©. At the centre of this halo 
. are two compact red massive galaxies. They are surrounded by hot gas, probably a 

l/^ ■ super- wind from merger-induced nuclear starbursts. This hot gas has shut down star 

CO \ formation in the non-nuclear region of these galaxies, leading to their red-and-dead 

Q.,^ . colours. 

■ A filament or lump of infalling cold gas is colliding with the hot gas phase and 
\ being shocked to high temperatures, while still around 30kpc from the red galaxies. 

The shock region is self-absorbed in Lyman-a but produces C IV emission. 

Further out still, the cold gas in a number of sub-halos is being lit up, most likely by 
a combination of tidally triggered star formation, bow-shocks as they plough through 

■ the hot halo medium, resonant scattering of Lyman-a from the filament collision, 
^ ' and tidal stripping of gas which enhances the Lyman-a escape fraction. The observed 

Lyman-a emission from the Blob is dominated by the sum of the emission from these 
sub-halos. 

On statistical grounds, we argue that Lyman-a blobs are not greatly elongated in 
shape, and that most are not powered by ionisation or scattering from a central active 
galactic nucleus or starburst. 

Key words: galaxies: haloes - galaxies: high-redshift. 



1 INTRODUCTION 

Sin ce their first discovery ([Francis et al ] |l996l : ISteideT et al. 
I2OO Q). Lyman-a blobs have held out the prospect of allowing 
us to study the spatially resolved flows of gas around mas- 
sive young galaxies. While some Lyman- a blobs surround 
high-redshift radio galaxies, and are clearly associated with 
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AGN jets (e.g. 'Reuland et al.' '2003', '2007ll. many are not. 
For these, a wide range of models have been proposed to ex- 
plain their extended morph ology and immense lu minosities, 
including superwinds (e.g. "Taniguchi et al."200l'), infalling 
gas (e.g. Francis et al. 1996; Dijkstra et al. 2006. ) and reso- 
nant scattering fe.g. ISteidel et al.ll201ll V 

Integral field spectroscopy is an ideal way to map the 
velocity structure in these low-surface-brightness objects. 
Such spectroscopy has been published for four Lyman- 
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a blobs to date, and all show a simila r pattern. LABI 
(|Bower et al ] |200^; IWeiimans etHI I2OIOI ) shows no large- 
scale systematic velocity field, and multiple velocity com- 
ponents. Many, but not all of the Lyman-a components 
seen here are seemingly associated with galaxies. lYang et aU 
(2011) observed blobs CDFS-LABl and CDFS-LAB2. They 
too find that the blobs break up into multiple velocity com- 
ponents, which seem to be associated with galaxies. They 
also find (for CDFS-LAB2) that the Lyman-a emission is 
self-absorbed at the systemic re dshift, with mo s t em ission 
coming at redder wavelengths. IWilman et all ^2005) ob- 
served LAB2 and saw absorption at one wavelength cross- 
ing their entire field, which they interpreted as absorption 
in a shell of gas swept up by a superwind, but an alter- 
native interpretation is that it is due to a self-absorbing 
clo ud of gas ionised from within, along the lines suggested 
bv lZheng &: Miralda-Escudd (|2002l ). Otherwise LAB2 shows 
the same multiple emission components seen in LABI. 

Thus all four of the Lyman- a blobs observed to date 
with spatially resolved spectroscopy show the same pattern: 
spatially extended but spectrally narrow emission compo- 
nents, often associated with stars, with no clear systematic 
velocity structure. 

In this paper we present integral field spectroscopy of 
the first Lyman- a blob foun d, (LAB 1_J2 1434423, hereafter 
Bl), originally discovered bv I Francis et al.l (|l996h . Our ob- 
servations have a spectral resolution at the wavelength of 
Lyman-a four times better than that of the previous obser- 
vations, which turns out to be unexpectedly informative. 

In an appendix, we also present spectroscopy of other 
Lyman-a blob candidates in the same region of the sky. 

Throughout this paper we adopt a cosmology with Hq = 
71 km s"^Mpc"\ Qm = .27 and Qa = 0.73, b ased on the 
five year WMAP results (|Hinshaw et al.ll2009l ). Distances 
are proper unless otherwise specified. 



2 TARGET 

The Lyman - g blob Bl, originally di scovered by 

I Francis et al. ([l996ll h a s also b een studied in [Francis et al.l 
(I1997II. [FYancis et al.1 (|200lh . IColbert et aU (|2006l ) and 
IColbert et al.1 (|201lh . This z = 2.38 blob remains one of the 
brightest and best-studied of the Lyman-a blobs. In this 
section, existing observations of this blob, as reported in 
the above papers, are summarised. 

Hubble Space Telescope (HST) imaging ([Francis et al.l 

[200l[ ) shows that Bl consists of two compact red galaxies 
separated by 0.8^^ (6 kpc), with a faint blue filament extend- 
ing around 50 kpc to the south (Fig[T|). The blue filament 
can be explained as a region of ongoing star formation, with 
little dust and a st ar- format ion rate of around lOMoyr"^ 
([Francis et al.[[2001). The two com pact red galaxies are more 



complicated. Francis et al.[ ([200l[ ) showed that their colours 
could be explained either by an old (> 750 Myr) stellar 

population, or by a younger one obscured by dust. 

The mid-IR photometry presented by IColbert et al.[ 

([201 l[ l allows this degeneracy to be resolved (Fig [2]). We 
combined the optimally weighted photometry , summed over 
all the compone nts, from [Francis et al.[ ([ 19961 ) , [Francis et al.[ 



|l997[) a nd Francis et al.l2001 ) w ith the Spitzer photome- 
try from IColbert et al.[ (|2006[ l and [Colbert et al.[ ([201 l[ l. to 



get the spectral energy distribution shown in Fig [21 As the 
different components of Bl were not resolved in most of 
these data, the spectral energy distribution is an integrated 
one. These data were then compared to spectral synthesis 
models from [Bruzual Charlot[ ([2003[ ). combined with dust 
absorptio n follow ing the empirical starburst extinction curve 
of Calzetti et al.[ (|l994). The compact red galaxies can be fit 
by a simple instantaneous burst stellar population of age of 
100-500 Myr, dust extinction with E{B-V) - 0.3 mag, and 
stellar mass of ^ 2 x 10^^ Mq. They are thus massive but 
"red and dead", with no substantive star formation for at 
least 10^ years. Additional blue emission from the filament 
explains the two shortest wavelength data points. Bl is a 
much more luminous source at 24 microns than this model 
predicts and shows poly-aromatic hydrocarbon (PAH) emis- 
sion. This implies that a dusty starburst with a star forma- 
tion rate of around 42OM0yr~^ is pres ent, most likely in th e 
core of one or both of the red galaxies ([Colbert et al.[[2006[ ) . 

The Lyman-a emission in Bl is genuinely diffuse. Inter- 
mediate band HST imaging shows that 13% of the Lyman-a 
emission comes from the eastern red compact galaxy, but the 
rest is not coming from clumps smaller than a few kpc in 
size, as otherwise we would have seen the clumps. 

Strong C IV (1 54.9 nm) emission is detected from Bl 
(I Francis et al ] [1996 1. Curiously, it is coming not from the 
location of the red galaxies, but from 25 kpc south along 
the filament, and occurs at a wavelength to the blue of most 
of the Lyman-a emission (Fig [3]). 

Bl lies at redshift 2.38, in the middle of a superclus- 
ter ([Francis Hewett[ [l993: Palun as et al.[[2004l l traced by 
QSO absorption lines and Lyman-a emitting galaxies. Spec- 
troscopy of three previously unclassified objects within this 
supercluster is reported in the Appendix A. 



3 OBSERVATIONS AND REDUCTION 

Bl was observed on the night of August 1 2011, using the 
WiFeS integral field spectrograph ([Popita et al.[[2007[ . [2oTol l 
on the 2.3m telescope at Siding Spring Observatory. WiFeS 
uses an array of 25 slits, each 38 arcsec long and 1 arcsec 
wide to obtain integral field spectroscopy of a 25 x 38 arc- 
sec field. The total integration time on Bl was 4.5 hours. 
Seeing was around 1.3 arcsec. The B3000 grating was used, 
giving a spectral resolution of 100 km/s. The data presented 
here were reduced using a custom set of scripts written in 
Mathematica. After bias and overscan subtraction, and fiat 
fielding, the individual slits were traced and straightened. 
Cosmic rays were identified by comparison with a median of 
the individual images. Very intense cosmic ray hits (proba- 
bly the result of beta-decay in the substrate) were found to 
produce an exponential tail, presumably due to some charge 
transfer problem. This appeared to be a threshold effect, 
only occurring for the brightest cosmic rays, and having a 
constant amplitude and shape for them. These tails were 
thus easily identified and removed. 

Wavelength calibration was done globally, to avoid in- 
troducing any row-to-row noise into the reconstructed data 
cube. Firstly, a polynomial wavelength dispersion relation 
was found for each spectral row in each slit. A fifth order 
polynomial was required to fit the arc lines in an individual 
row. Each of the derived polynomial coefficients was then 
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Figure 1. The left panel shows Hubble Space Telescope F450W image of Bl (greyscale) overlaid with the integral Lya flux measured 
in this paper (contours, smoothed by a Gaussian with standard deviation 0.5^^). Note that previous observations showed that the Ly-a 
emission continues at low surface brightnesses much further out than shown here. The right panel shows a sketch of our model for Bl. 



plotted as a function of slit number and position within each 
slit. A linear or cubic fit to the variation in each polynomial 
coefficient was then derived. These fits were then used to re- 
construct the wavelength solution for any part of the data. A 
root-mean-squared residual of O.OlSnm was achieved across 
all arc-lines in all parts of the data, with no systematic pat- 
terns to the residuals. Sky lines in reconstructed data cubes 
were well aligned both along and across slits. 



4 RESULTS 

The only line detected in the WiFeS data was the Lyman-a 
line of blob Bl. Different views of the data cube centred on 
Bl's Lyman-a line are shown in Figs (H O [6] and O 

The velocity field is clearly a complex one. A number 
of conclusions can be drawn: 

(i) The Lyman-a emission from Bl is elongated by around 
80 kpc in a north-south direction. The northern end corre- 
sponds to the red galaxies, and the southern end to the blue 
filament. 

(ii) The peak emission varies as a function of wavelength 
between the northern and southern ends, giving it a zig-zag 
appearance in velocity space (FigsjSlHj). 

(iii) There is a gap in the emission at around 411.7 nm 
wavelength. 

(iv) The emission, particularly at the southern end, shows 



a complex velocity structure, with multiple emission and/or 
absorption components which are spectrally very narrow (< 

100km s"^). 

(v) The C IV emission seen bv lFrancis et al.l(|l996l l comes 
from the southern end of Bl, but at a wavelength blueward 
of nearly all the Lyman-a emission. 

Note that these observations are only picking up the 
higher surface-brightness parts of Bl: our previous observa- 
tions showed that low surface-brightness Lyman-a emission 
extends far to the north-east of the region discussed here. 
We also know that a QSO sight-line passing around 160 
kpc to the north-west intercepts a Ly man- limit absorption- 
line system ([Francis k Willigerll2QQ4 ). indicating that dense 
neutral hydrogen can be found far beyond the extent of the 
observable Lyman-a emission. 

This zig-zag velocity structure, interrupted by a gap, 
cannot easily be explained by a single rotational, infall 
or outflow model. It could be explained by a combina- 
tion of separate components. Note that it is unlikely that 
the velocity structures we see are related one-to-one with 
gas structures. Due to the high optical depth in Lyman- 
alpha, the emission from any slab of gas structure is 
likely offset in velocity from its systemic velocity (e.g. 
IZheng k Miralda-Escudl l2002l : ISchaerer etaD l201lh . and 
one slab of gas may produce two or more emission compo- 
nents, typically in the velocity wings. Nonetheless, radiative 
transfer effects typically broaden lines, so our upper limit 
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Figure 3. Slice through the data-cube, of width one arcsecond. Wavelength is shown horizontally and position north/south vertically. 
Spatial pixels are one arcsec wide, spectral pixels are 0.075 nm wide. The data have been spatially smoothed as in Fig[5l but not spectrally 
smoothed. 



on the velocity dispersions of around 100 km s~^ is thus an 
upper limit on the velocity dispersion in the gas slabs from 
which this emission arises. 



with galaxies. However, our better velocity resolution, com- 
bined with the extensive multi- wavelength data-set available 
for this well studied blob, allows us to draw some somewhat 
more detailed conclusions than previous authors. 



5 DISCUSSION 

Our observations show a very similar pattern to the previ- 
ous integral field observations: no simple large-scale system- 
atic velocity field, multiple sub-components, some associated 



Our discussion starts with a puzzle, an hypothesis and 
a constraint. 
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Figure 4. Three dimensional view of the data cube, with the region of strongest flux shown. 



5.1 Puzzle 

The puzzle concerns gas in Bl. The extensive Lyman-a neb- 
ula implies the presence of large amounts of cool gas. In con- 
trast, the red-and-dead colours of the compact red galaxies 
imply that there has been little star formation over a period 
lasting at least 100 Myr, and hence little cool gas. On the 
other hand, the mid-IR luminosity and PAH detection in- 
dicate prolific star formation and hence large quantities of 
molecular gas. 

The resolution to this puzzle is to have different phases 
of gas on different scales. The starburst is most likely very 
compact (< 1 kpc) and found in the nuclear regions of one 
or both of the compact red galaxies. It is so obscured by 
dust that no sign of it is seen at rest-frame wavelengths 
shorter than around 2.5/im (Fig [2]). This nuclear starburst 
(or starbursts) is probably driven by the interaction between 
the two red galaxies, which lie at a projected separation of 
only 6 kpc. This pict u re is co nsistent with the expectations 
of the iKennicutt Jrl (|l998bh star formation law in which 
the surface rate of star formation is related to the surface 
density of gas though Ssfr oc Xlgas- Regions of high specific 
star formation are therefore regions of high gas (and dust) 
surface densities, and are therefore likely to be heavily dust- 



obscured. The Lyman-a emission we detect from one of the 
red galaxies would presumably be coming from the outskirts 
of the galaxy. 

On scales of a few kpc, the gas is probably too hot to 
allow star formation. This explains the red-and-dead colours 
of the compact red galaxies. A superwind is likely present, as 
nuclear starbursts as luminous as that in Bl v irtually always 
produce superwinds in the local universe fe.g. Veille ux et aP 
.2005 ) . This could have swept out any cool gas and trun- 
cated star formation i n the outer regions of these galaxies. 
^Taniguchi et al.l (|200lh suggested such superwinds as power- 
sources for Lyman-a blobs. The age of the stellar population 
('-^ 10^ years) is comparable to the merger timescale, and the 
current star- format ion rate, acting over 10^ years would pro- 
duce a stellar mass comparable to that observed. The gas 
may also, however, h ave been he ated by a concealed active 
galactic nucleus (e.g. iKing et aD 2011) or by the liberation 
of po tential energy from infalling gas (e.g. iJohansson et all 

We now turn to the neutral hydrogen responsible for 
the Lyman-a nebula. This presumably lies further out, and 
is the subject of our hypothesis. 
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Figure 5. Images of Lyman- alpha emission from Bl, as a function of observed- frame wavelength. Each pixel in these images is one 
arc-second square. North is to the top and East to the left. The images have been smoothed by a Gaussian with a standard deviation 
of one pixel (one arcsecond). Each slice covers a wavelength range of 0.075 nm. The box in the top-left image shows the field of view of 
Fig El 



5.2 Hypothesis 

Our data show that the extended Lyman-a flux is typically 
tens of kpc from the red galaxies (which we assume lie at 
the centre of the dark matter halo), spatially extended on 
scales of several kpc, but breaks up spectrally into narrow 
sub-components. All these properties sound similar to those 
of a family of dwarf galaxies surrounding a massive galaxy. 
Which led us to the following hypothesis: that much of the 
extended Lyman-a flux is coming from neutral gas within 
dark-matter sub-halos, moving within the much more mas- 
sive dark matter halo centred on the red galaxies. 

How feasible is this model? The relative velocities of 
the sub-clumps are ^ 700km s~^, and they are spread over 
^ 50kpc, which gives a virial mass of the main halo of 
^ 10^^ M0. We selected ten s uch massive halos from the M il- 
lenium II n-body simulation (|Bovlan-Kolchin et al.ll2009l l at 
this redshift, and looked at the sub-halo population of these 
ten big halos. There were typically 10 - 20 sub-halos within 
each big halo. They were spread over around 100 kpc, and 
had a line-of-sight velocity dispersion of around 500km s~^. 
And the velocity dispersions of the individual sub-halos were 
typically around 70km s~^. Thus, the observational and the 
predicted parameters appear to match each other very well. 
These sub-halos in the simulation had typical dark matter 



masses of ^ 10^^ or more and half-mass radii of 1 - 5 
kpc. 

Would one expect sub-halos such as these to contain 
neutral hydrogen? We picked out a sample of small halos 
with masses of - 10^^ Mo from the semi-analytic models of 
IPe Lu cia & Blaizot (2007) applied to the Millenium I simu- 
lation (jSpringel et al.ll2005l ) , and they typically had cold gas 
masses of - 10^ M©. So current semi-analytic models do sug- 
gest that at these redshifts, halos with masses of — 10^^ M© 
will contain reasonable amounts of neutral gas. 

Why are these sub-halos radiating Lyman-a emission? 
The sub-clumps we observe have Lyman-a luminosities of 
— 2 X 10"^^ VK, which are well within the range found in ob- 
served Lyman-a emitting galaxies. However, the space den- 
sity of — 10^^ M© halos is vastly greater th an the space 
densi ty of Lyman-a emitting galaxies (e.g. iGuaita et al.l 
I2OIOI I. Thus most such halos do not radiate at this level. 
The equivalent widths are also different: we compared the 
equivalent widths of blobs with those of compact Lym an- a 
emitting galaxi e s in t wo surveys: IPalunas et aP (|2004l l and 
iMatsuda et al ] (|2004l l. Both surveys measured the equiva- 
lent widths (or equivalently narrow-band minus broad-band 
colour) for both the blobs and other Lyman-a emitting 
galaxies in their fields, and in both cases, blobs were picked 
purely on the basis of size in the narrow-band image. In 
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Figure 6. Grid of spectra for individual spatial (1x1 arcsec) pixels. North is to the top and East to the left. The data were spatially 
smoothed as described in Fig [5] The location of this grid is shown in FigO 



both cases, Lyman-a blobs tend to have higher equivalent 
widths. A Kolmogorov-Smirnov two-sample test shows that 
the difference is significant at the 95% confidence level for 
the Palunas data, and at the 99.7% confidence level for 
the larger Matsuda dataset (using objects brighter than a 
narrow-band magnitude of 24.2). We therefore conclude that 
Lyman-a blobs have significantly higher equivalent widths. 

Thus for the hypothesis to be valid, we need some mech- 
anism which increases the Lyman-a luminosity and equiva- 
lent width of ^ 10^^ Mq sub-halos when they are within the 
dense environment of massive halo, compared to the values 
seen in similar mass halos in less dense environments. We 
will discuss the possibilities in Section [5.41 

Note that some of these sub-halos may contain stars 
at this redshift. Sensitive observations may thus find that 
most, but probably not all of th e Lyman-a cl umps have 
faint galaxies within them. Prescot t et al.l (|201 1a) recently 
presented imaging showing that one Lyman-a blob does in- 
deed seem to comprise a compact group of dwarf galaxies, 
though in this case the galaxi es were not clea r ly ass ociated 
with the Lyman-a emission. Weijmans, et all ((2OIQI ). how- 
ever, did find that most clumps of Lyman-a emission were 
centred on faint galaxies. So this is plausible. 



5.3 The C IV Constraint 

We now come to our constraint: the strong C IV (154.9 
nm) emission detected comi ng from ^ 25 kpc south of the 
red galaxies. As discussed in iFrancis et al.l (|l996 ), C IV re- 
quires a hard ionising spectrum and must therefore either 
come from Active Galactic Nucleus (AGN) photoionisation 
or from fast shocks. We will discuss these possibilities in 
turn. 



5.3.1 AGN Model 

For the AGN model to work, we would require that at least 
one of the red galaxies hosts a type-2 AGN. Our line-of-sight 
to the AGN is blocked by a dusty torus, but UV radiation 
is escaping at least in the southerly direction, and impact- 
ing on some relatively dense gas cloud, w here it produces 
Lyman-a and C IV via photoionisation fe.g. Haiman ReesI 
2001; Adelberger et al. 2006; Francis & McDonnell 2006). A 
second mechanism, excitation by collision with a radio jet, 
i s ruled out by our upper limit on the radio flux from Bl 
(jFrancis et al.lll996l ). 

The first argument against this model is the detection 
of PAH emission from Bl, which suggests that if a con- 
cealed AGN is present, its energy output is dominated by 
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that from star formation. Note however that the PAH emis- 
sion was only tentatively detected (Col bert et ali ,2Qll), and 
that even an AGN which makes no significant contribution 
to the total energy output could still produce significant 
photoionisation. 

A second argument against the AGN model is statisti- 
cal. If the AGN model is correct, and the ratio of obscured 
to unobscured AGN is r, then if we see n blobs in a field, 
we should see n/r unobscured AGN. We know that r is 
not significantly n e ater than 3 (e.g. [ Webster et aD Il995l : 
iFrancis et alll2004al : [Xreister et al.ll2009f L 

In the Palunas et al. field, we have three (or four if 
you count B5) Lyman-a blobs. How many QSOs lie in this 
field at our redshift? Two QSOs were identified at the clus- 
ter redshift, one wit h B — 21.24 a t coordinates 21:42:31.33 
-44:30:16.8 found bv lFrancis et al.l (|2004ah . and B29 identi- 
fied by our GMOS observations described i n Appendix A, 
with B — 23.16. In the field surveyed by iMatsuda et al.l 
there are 35 Lyman-a blobs, and five known QSOs 
with optical magnitudes in the range 20.5 to 24.5, found 
in the NA S A Ex tragalactic Database. In the survey by 
I Yang et"aL 1 iooi), there are four Lyman- a blobs and no 
QSOs down to the Sloan Digital Sky Survey limit of ^ < 
19.4. 

Combining these samples, we find that we can get suit- 
able numbers of unobscured QSOs to be consistent with the 
blob population, but only if we include faint QSOs: those 
down to y ~ 23 or fainter. 

Could QSOs this faint feasibly enhance the Lyman- 
a emission from nearby gas clouds to give the luminosity 
we observe in Lyman-a blobs? We can answer this em- 
pirically by looking at the observed ratio of QSO mag- 
ni tude to the flux of extended Lyman-g emission around 
it. IWeidinger et al.l (|2005l ) and IChristensen et al.l (j2006h be- 
tween them detected extended Lyman-a emission around 
9 QSOs. The average extended Lyman-a flux was 3.2 x 
10~^^erg cm~^s~^, which is comparable to Lyman-a blob 
fluxes, so QSOs certainly can produce extended emission 
as seen in Blobs. However, the QSOs were very bright: av- 
eraging V — 18.2: much brighter than the QSOs found in 
the various Blob fields. And many QSOs this bright in the 
Weidinger and Christensen samples had no detectable ex- 
tended Lyman-a flux around them, so this average ratio of 
extended flux to nuclear flux is really only an upper limit. 
If we assume that the extended flux scales with the nuclear 
flux, then the QSOs in the blob fields are around two orders 
of magnitude less luminous than would be needed to pro- 
duced the observed Lyman-a blob flux. Would we expect 
the extended flux to be proportional to the photoionising 
luminosity? It is possible that the most luminous QSOs ac- 
tually suppress the Lyman-g emission o f nearby galaxies, 
as suggested by iFrancis Bland- HawthornI (|2004l l . In this 
case, less luminous QSOs might be surrounded by stronger 
Lyman-a fuzz than the most luminous ones (which are the 
only ones studied to date), and hence might be plausible 
unobscured counterparts of Lyman-a blobs. We note, how- 
ever, that we did not detect extended Lyman-a emission 
from either of the QSOs lying within our supercluster. 

We therefore conclude that if blobs contain concealed 
QSOs, and these concealed QSOs are of comparable lumi- 
nosity to the other QSOs seen in the blob fields, then they 
are unlikely to be luminous enough to produce the extended 



Lyman-a emission through photoionisation. If the concealed 
QSOs are much more luminous, then it becomes difficult to 
understand why we don't see at least a few comparably lumi- 
nous unobscured QSOs in these fields, unless the radiation 
from these QSOs is completely blocked, in which case it can- 
not escape to photoionise the extended Lyman-a emission. 
Thus while the AGN model cannot be ruled out for indi- 
vidual blobs, for statistical reasons it is unlikely that AGN 
photoionisation is the dominant power source of the popu- 
lation of Lyman-a blobs. 

5.3.2 Fast Shocks Model 

IFrancis et al.l (| 19961 ) and IFrancis et al.l (|200lh originally sug- 
gested that fast shocks could cause the C IV emission from 
Bl, and demonstrated that they could produce ample line 
emission provided the density of the medium into which they 
were passing was at least 10~^cm~^ and the shocks had be- 
come radiative, which required a timescale of ^ 10^ years, 
which is comparable to the infall time for Bl and hence 
plausible. 

Fast shocks are almost inevitable given the high velocity 
dispersion and observed lumpiness of the gas around Bl. A 
plausible location for the shocks would be a collision between 
a superwind from the red galaxies and gas infalling along a 
filament from the south. 

What is the infalling gas which is being shocked? One 
possibility is gas within an infalling sub-halo, perhaps an 
unusually massive one. Another is cold-mode infall (e.g. 
'Faucher-Gigue re et al.l l2010l : ICoerdt et aP l2010l ). The in- 
falling gas cannot, however, be primordial, as we need it to 
contain carbon. The presence of a superwind is not neces- 
sary for this model to work: the impact of an infalling cold 
gas cloud on the hot central medium will suffice to drive 
suitable shocks. 

5.4 Making the Sub-Halos Glow 

Our model, based on the puzzle, hypothesis and constraint 
discussed above, is summarised in Fig[T] To make it work, 
however, we have to find some mechanism to make the gas 
in each sub-halo emit a Lyman- a luminosity of > 2 x 10^^ VK. 
There are several plausible mechanisms and we will discuss 
them in turn. 

5.4' i Induced Star Formation 

The ram pressure of their passage through the hot medium, 
or tidal interactions, could perhaps trigger star formation in 
these sub-halos, which would in turn generate the Lyman-a 
emission. If we take the observed Lyman-a luminosities, as- 
sume a typical unobscured ratio of Lyman-a/Ha of ^ 8 : 1, 
and apply the relati on between Hq lumi nosity and star 
formation rate from iKennicutt Jrl (|l998al ). then we infer 
star formation rates in each sub-clump of around 2-3 so- 
lar masses per year, and hence multiplying by the number 
of sub-clumps, a total star formation rate of around 10- 
20 solar masses per year. This would just about be con- 
sistent with producing the observed UV emission from the 
blue filament (which gave an inferred star formation rate of 
around lOMoyr"^). In practice, however, the optical depth 
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in Lyman-a is likely to be so high that ratios of Lyman- 
a/Ha will be closer to ^ 1 : 1, which would raise the in- 
ferred star formation rate to levels that are inconsistent with 
the faint ness of the observed blue continuum flux. Another 
way to say this is that the equivalent width of the Lyman- 
a emission is higher than that observed for other Lyman-a 
emitting galaxies, as discussed in Section 15.21 Star forma- 
tion is thus unlikely to produce all the Lyman-a emission, 
though it could produce some of it. 

5.4-2 Inf all-powered emission 

A typical sub-halo, according to semi-analytic models, will 
have a cold gas mass of > 10^ M©, and an infall time of 
^10^ years. If we divide the gravitational potential energy 
of the gas in a sub-halo by the infall time, we get a power 
of ~ 10^^ VK, which is 500 times greater than the observed 
Lyman-a luminosity. Thus if we can convert even 1% of the 
gravitational potential energy into radiation, we can explain 
the brightness of the sub-halos. 

One mechanism for doing this (e.g. IPiikstra Loebl 
I2OO 9) is cooling radiation of thermally heated gas. Another 
is shock ionisation. If we assume that each sub-halo has a 
neutral gas surface area of several tens of square kpc, then 
a Lyman-a luminosity density of a few xlO~^Wm~^ is 
needed to reproduce the observed luminosities. If a shock- 
wave is driven into this surface as the sub-halo moves 
through the ambient gas, this shock will prod uce a diffuse 
Lyma n-a emission. The model shock grids of en et al. 
(2OO&) show that such a Lyman-a luminosity density can be 
easily produced, provided pre-shock densities are > O.lcm"^ 
and shock velocities are > 200km s~^. Higher velocities or 
densities give much higher luminosities, allowing smaller gas 
surface areas, or more self-absorption. 

Thus infall is a viable model for producing at least a 
substantial part of the observed Lyman-a luminosity. 

5.4-3 Scattered Lyman-a from a central source 

One possibility is that the emission we see from the sub-halos 
is actually produced elsewhere, and is just being resonantly 
scattered by the sub-halos. The emission could be coming 
from the same region that generates the C IV emission (ac- 
cording to our model, a region ionised by fast shocks where 
infalling gas hits a superwind) , or escaping from the nuclear 
starbursts in the red galaxies. We do not see Lyman-a emis- 
sion directly from these regions, due to their high optical 
depths, but instead we see part of it scattered into the wings 
of the velocity distribution and off nearb y gas clouds. This 
mechanism was invoked by Steid el et al.l ([2OII) to explain 
the faint diffuse Lya halos seen in stacked images of high 
redshift g alaxies, and can lead t o polarised Lyman-a emis- 
sion fe.g. iDiikstra Loebl lioosl ). iHaves et all (|201lh found 
that th e Lyman-a em i ssion in LABI was indeed polarised, 
though [ Prescott et al.1 (|2011bl) did not detect such polarisa- 
tion in much shallower observations of a different Lyman-a 
blob. 

How plausible is such a model in the case of Bl? 
ISteidel et al. (2011) found that the extended Lyman-a flux 
typically exceeded the compact central line flux by a factor 
of ^ 5, and in our case, the extended Lyman-a flux exceeds 



the central Lyman-a flux (from one of the red galaxies) by 
a factor of 7, which is comparable. 

If scattered Lyman-a emission originally came from the 
superwind/infalling gas collision, then the UV emission pro- 
ducing it came from fast shocks (as we know from our C IV 
detection) and would hence have a very hard spectrum. This 
means that we would expect very high equivalent widths and 
little detectable UV continuum, as observed. 

If, on the other hand, the scattered Lyman-a emis- 
sion originally cani e from star formation, as proposed by 
ISteidel et all (|201lh . we might expect to see UV emission 
from this star formation. We can est imate how b r ight it 
should be by scaling from the results of lSteidel et al ] (|201lh . 
Bl has a Lyman-a surface brightness in its central few arc- 
sec which is ~ 30 times that of the halos seen by Steidel 
et al., so one would expect blue continuum emission ^ 3.5 
magnitudes greater than the Vab ^ 25 of the galaxies stud- 
ied by Steidel - i.e. we should expect V ^ 21.5. In fact, 
Bl has V ^ 23.5. The red galaxies in Bl show no evidence 
in their optical colours of any ongoing star formation: we 
only see evidence for it in the mid-IR. A similar problem 
occu rs in LABI, where th e probable central galaxy is also 
red (IChapman et al.ll2004l ). In both Bl and LABI there is 
plenty of star formation taking place, which is seen in the 
mid-IR and sub-mm, but little or no detectable UV light. 

To put this another way, the star formation rate in- 
ferred from the observed blue light in Bl (the blue filament) 
of ^ lOM0yr~^ is insufficient to produce the phenomenal 
observed Lyman-a luminosity. We would need to tap into 
the much greater star formation rate which is concealed by 
dust, most likely in the red galaxies (^ 42OM0yr~^). 

The only way out of this would be to make the radiation 
from starbursts anisotropic. We would not see it from our 
orientation, but it would escape in other directions, analo- 
gous to unified models of AGN. If this was the case, we would 
sometimes expect to see blobs in which the central starburst 
was exposed. These would look like compact Lyman-a emit- 
ters surrounded by relatively faint fuzz, rather than diffuse 
blobs, but should have Lyman-a luminosities greater than 
that of blobs. There are two potential exposed starbursts in 
the Palunas et al. su percluster: B5 (Appendix A) and B40 
([Francis et al.l l2004bY Both have Lyman-a luminosities and 
line widths comparable to the three blobs in this field. B5 has 
a compact nucleus surrounded by relatively faint fuzz, while 
B40 is spatially unresolved. But the extended flux of these 
objects is much fainter than that of the Lyman-a blobs, so 
if viewed from some other orientation they might well be 
blobs, but much fainter ones. There are thus no possible 
counterparts of the most luminous blobs in this field. 

Thus the scattering model is feasible, but only if we are 
scattering emission generated by fast shocks, rather than by 
star formation. 



5.5 Tidal streams 

Regardless of the mechanism or mechanisms producing the 
Lyman-a emission from the sub-clumps, one would expect 
some tidal stripping of gas from them in this dense en- 
vironment, perhaps producing Magellenic- stream style gas 
streams. This extended narrow gas morphology will increase 
the escape fraction of Lyman-a and help enhance the ob- 
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served brightness of these sub-halos regardless of the cause 
of the emission. 



5.6 Intrinsic shape of the Sub-halo Swarm 

If Lyman-a blobs are collections of sub-halos, are these 
collections roughly spherical or are they elongated? There 
i s some ev idence for blobs lying in filamentary structures 
(|Erb et al.lE)lli) but is this reflected in their shapes? 

In this section, we develop a simple toy-model of elon- 
gated Lyman- a blobs, taking into account a selection effect 
which biases samples towards finding end-on blobs. We com- 
pare this to d ata on the apparent shapes of Lyman-a blobs 
from iMatsuda et al.l ( 2004), and show that blobs are not 
highly elongated on average. 

Let us model Lyman-a blobs as cigar-shaped ellipsoidal 
distributions of line-emitting clouds. To keep the model sim- 
ple, we assume that the observed Lyman-a surface bright- 
ness of a blob is proportional to the length of our line of sight 
through the ellipsoidal distribution. We further assume that 
the long axis of the blob is a factor r longer than the mi- 
nor axis, and that we are observing at an angle to the 
long axis. Simple geometry then tells us that the observed 
surface brightness b is given by: 



6oc V'l + (r2-l)cos2(6'). 



(1) 



Any narrow-band search for Lyman-a blobs will have a sur- 
face brightness cut-off. Blobs that would normally lie below 
this cut-off will make it into the sample if viewed end-on, so 
any sample will be biased towards blobs seen at small values 
of 0. How big is this bias? To estimate this, we need to know 
how the space density of Lyman- a blobs increases as the 
surface brightness decreases. To parameteris e this, we plot- 
ted t he cumulative number of blobs found bv lMatsuda et al.] 
(^4) against their surface brightness, down to the rough 
completeness limit of this sample (Figl?)). We fit the cumula- 
tive space density in this survey with a straight-line model: 
the best fit is: 

log^Q (Cumulative Number of Blobs) oc 1.0 mag^^ (2) 

Thus if blobs seen end-on are a magnitude brighter than 
those seen side-on, the end-on ones will be over-represented 
in a sample by a factor of around 10. 

We assume that the cumulative number density has a 
power-law dependence on limiting magnitude, with an in- 
dex chosen to reproduce this factor of ten. For any given 
surface-brightness enhancement, we can then calculate how 
many more blobs we should see. We further assume that 
in the absence of any selection effects, blobs are randomly 
orientated in three dimensions, which means that the prob- 
ability of seeing a blob from a given angle is proportional 
to sin(^). The observed distribution is then the product of 
the random orientation with the bias towards end-on blobs, 
which comes out as: 



sin 6* X 10 



(1.0x2.5 log 10^1 + (r2-l) cos2 6) 



(3) 



which simplifies to 



(l+ (r'-l)cos'6')'''sin6' (4) 

which is plotted in Fig[8l 

Thus even quite modest elongations will have a strong 



impact on the viewing angle at which Lyman-a blobs are 
seen. We thus confirm that there can be a bias towards small 
viewing angles 0. 

If we assume a power-law intrinsic distribution of elon- 
gations (r values) , we can then do a monte-carlo calculation 
of the distribution of observed axial ratios projected onto 
the sky (the ratio of the observed long- to short-axis of the 
blob) , and compare this wit h the observed distribution taken 
from iMatsuda et al.l (|2004l ). The results are shown in Fig[9l 

This figure shows that most of the Lyman- a blobs ob- 
served by Matsuda were not greatly elongated. The best fit 
is where the number of blobs with a given intrinsic elonga- 
tion [R) is proportional to . Models in which elongated 
blobs are more common over-predict the number of blobs 
with large observed elongations. 

The Millenium II simulations analysed in sect ion [5?2l are 
also not strongly elongated. Some of the massive galaxies lie 
at the intersection of multiple filaments, but because the 
sub-halos lie along more than one filament, the projected 
distribution is typically not greatly elongated. 

We therefore conclude that Lyman-a blobs are not in 
general strongly elongated (though clearly a few are). Ob- 
served samples are biased towards viewing elongated blobs 
down their long-axis, but this bias is not sufficient to explain 
the observed predominance of nearly circular blobs. 



6 CONCLUSIONS 

As our observations of Lyman-a blobs become better, 
the blobs appear more complicated. There is already ev- 
idence for multiple energy sources in different blobs (e.g. 
Colbert et al. 2011). This paper suggests that for one blob at 
least, the origin of the surrounding gas is also complicated, 
combining merger-driven dense molecular gas in the nuclei, 
hot gas perhaps from a superwind and infalling lumpy colder 
gas. Our observations are consistent with seeing a massive 
galaxy being assembled by minor mer gers , as des cribed by 
iKereg et al.l (|2005|), I Johansson et al.l ((2009|) and lOser et al.l 
(|201lh . though with a major merger taking place at the same 
time. 

Cosmological hydrodynamical models which include 
some combination of wind s and infall are rapidly increas- 
ing i n sophistication ( e.g. Zheng et al.l 20ld: ICoerdt et"al] 
'2010"; Schaere r et al.l 1201 ll : iFaucher-Giguere et al.l |2oTq[ 
Zheng et al. 20111) but most fail to predict the narrowness 
of the velocity sub-structures we see. The exception are the 
models of Go erdt et al. but the relative narrowness 

of the emission they predict is likely an artefact of the lack 
of radiative transfer in their models. 

We conclude that high spectral resolution observations 
are providing an interesting challenge to the models. Few 
instruments are set up to make these high resolution obser- 
vations in the blue, but those that can should pursue this. 
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Figure 7. The cumulative number of Blobs found by Matsuda et al, as a function of narrow-band magnitude. Each point represents one 
blob. The line is a best fit model. 
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Figure 8. Probability of observing a Lyman-o; blob as a function of orientation angle 0, for three different elongations r. 
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Figure 9. Comparison between the observed (histogram) and modelled (points) distributions of observed axial ratios. The two models 
are for populations of blobs with power-law distributions of r, with indices -1 and -3. 



APPENDIX A: SLIT MASK SPECTROSCOPY 
OF OTHER CANDIDATE BLOBS IN THE 
FIELD 

In this appendix, spectroscopy is presented of three other 
candidate Lyman-a emitters which lie in the same super- 
cluster as Bl. All lie several arcminutes from Bl and so do 
not directly impact on its physics. Their properties are, how- 
ever, used in the various statistical arguments in the main 
body of the paper. 

We had previously obtained spectra of most bright po- 
tent ial Lyman-o; emitte rs in the supercluster surrounding 
Bl (jFrancis et al.ll2004bl l. but three needed further observa- 
tions. In this section we present spectroscopy of these three. 

One, B5 in the list of Palunas et al. (2004) was classified 
as a potential Lyman-a blob, albeit an unusual one. It con- 
sisted of a luminous narrow-band emitting core surrounded 
by faint narrow-band excess fuzz. Our previous s pectrum 
of this object was ambiguous ([Francis et al.ll20"04bl l: we saw 
a potential Lyman-a line but also two other possible weak 
lines which did not m atch. The second object, identified by 
IPalunas et aD (|2004l ) and called B29 was the most lumi- 
nous compact potential Lyman-a emitter in the field, but 
had not been previously observed due to fibre collision con- 
straints with the spectrograph. It is located at 21:43:05.90 
-44:27:21.0 (J2000). The third, B30, was a fainter narrow- 
band excess object found by Palunas et al, also not observed 
previously due to fibre placement issues. It is located at 
21:43:06.42 -44:27:00.6 (J2000). 

B5, B29 and B30 were observed with the CMOS spec- 



trograph on the Gemini South Telesco pe. The obser v ationa l 
set-up with CMOS is described by Scar lata etHI (|2009h : 
both objects were observed simultaneously for an exposure 
time of 8,800 seconds, in 0.6 arcsecond seeing, on the night 
of 2005 October 9th. 

From B5, we detect a single very strong emission line 
at 410.0 nm, with a velocity width of 100 0km s~^. Th e ten- 
tative additional lines seen by F rancis etHI (|2004bl ) were 
ruled out. On the basis of the line width, extended emis- 
sion and equivalent width, so the line must be Lyman-a. 
The Lyman-a emission was slightly extended along the slit, 
so we confirm that B5 is an extremely luminous extended 
Lyman-a emitting object at 2; = 2.38. 

B29 turned out to be a QSO at the supercluster red- 
shift. It showed strong broad Lyman-a peaking at around 
412.0 nm, as well as a number of other typical QSO broad 
lines. The line widths were around 440 nm (full width at half 
maximum height). Strong associated absorption was seen in 
Lyman-a and C IV (154.9 nm), slightly red-shifted with re- 
spect to the emission-line peak. 

B30 showed a single strong emission line at 412.2 nm, 
with a velocity width of 900km s~^, and is hence confirmed 
as a Lyman-a emitting galaxy in the supercluster. 
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